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ABSTRACT 

We develop, implement and characterise an enhanced data reduction approach which deliv- 
ers precise, accurate, radial velocities from moderate resolution spectroscopy with the fibre-fed 
VLT/FLAMES-I-GIRAFFE facility. This facility, with appropriate care, delivers radial velocities 
adequate to resolve the intrinsic velocity dispersions of the very faint dSph dwarf galaxies. Impor- 
tantly, repeated measurements let us reliably calibrate our individual velocity errors (0.2 < 5v < 5 
km s~^) and directly detect stars with variable radial velocities. We show, by application to the 
Bootes I dwarf spheroidal, that the intrinsic velocity dispersion of this system is significantly below 
6.5km/s reported by previous studies. Our data favor a two-population model of Bootes I, consisting 
of a majority 'cold' stellar component, with velocity dispersion 2.4j^Q g km/s, and a minority 'hot' 
stellar component, with velocity dispersion ^ 9km/s, although we can not completely rule out a sin- 
gle component distribution with velocity dispersion 4.6*^0 6 '^^Z^- We speculate this complex velocity 
distribution actually reflects the distribution of velocity anisotropy in Bootes I, which is a measure of 
its formation processes. 

Subject headings: methods: data analysis — techniques: radial velocities — galaxies: dwarf — galaxies: 
individual (Bootes I) — galaxies: kinematics and dynamics 



1. INTRODUCTION 

There is continuing interest in analysis of the num- 
ber, nature, masses, and evolutionary histories of the 
dwarf spheroidal (dSph) satellite galaxies, found in mod- 
erate numbers around both the Milky Way Galaxy and 
M31 in the Local Group. They have typical half-light 
radii greater than lOOpc, low surface brightnesses (^ 
25-30 mag/sq. arcsec), and central velocity dispersions 
of several km/s, implying that the luminous compo- 
nent is embedded in a dominant extend ed dark mat- 
ter halo (M/Ly ^ 10 - 1OOM/Lv,0; (jMateol [19981: 
iGilmore et al.l l2007f) ) . Their chemical abundances are 
low in the mean, correlate with dSph system lumi- 
nosity, show real intrinsic dispersion, and have chem- 
ical element ratios systematically different from those 
of Galactic halo field stars over the metallicity range 
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covering most Galactic halo stars (iKirbv et al.l [20091 : 
iNorris. Wvse. Gilmore eraI|[20Tol: iTolstov et arn2009[) . 
Their stellar populations differ systematically from those 
in the Galactic field. All the astrophysical evidence 
shows they are the oldest surviving bound systems, prob- 
ably forming very early frorn purely primordial gas, e.g., 
INorris. Wvse. Gilmore et al.l poiol ). How do they relate 
to the very large numbers of surviving dark-matter halos 
predicted by standard structure formation models? 

The most luminous dSphs around the Milky Way, 
the 'classical' dSphs, were discovered photometrically 
through the second half of the twentieth century, apart 
from the nearest and largest, the Sgr dSph, which 
was di scovered in position-yelocit y -photo metry phase 
space ([Ibata. Gilmore. IrwinI 11994 |1995[) . More re- 
cently, following availability of the large-area pho- 
tometric data fro m the Sloan Digital Sky Survey 

(SDSS; lYork et ahl '2000') . a three-times larger sample 
of dSph has been discovered, primarily of dSphs with 
lower intrins ic luminosities, extending to the 'ultra- 
faint' dSphs (iWillnian et al.|[2005l: iBelokuroy et al|[2006l: 
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20071 iBelokurov et all I2008L I2009D These objects 
have extremely low surface brightnesses (down to 30 
mag/arcsec^) and low luminosities (down to Ly ^ 
lO^Ly_0), such that the presence of a single giant star can 
substantially affect the luminosity of the entire galaxy 
([Martin et al.l [200g). 

Knowledge of the numbers and masses of dwarf satel- 
lites remains crucial for understanding local galaxy for- 
mation, since high-resolution cosmological N-body sim- 
ulations of galaxy formation generally predict a num- 
ber of dark matter 'sub-haloes' that is still an order 
of magnitude larger than the number counted in and 
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extrapola ted from observations (the "Missing S atellites 
Problem" (|KlvDin et alJUOOl IMoore et alJll999D '). Con- 
siderable efforts have been expended i n simu l ating dwarf 
galax y formation (Barkana & L oeb| | 1999t iSomerviH^ 
I2002t iBenson et aT1l2003t IRkottl |20100, attempting to 
lessen the tension bet ween predictions from simulations 
and the observations (iKoposov et alll2009l : iMaccio et alj 
[2OT0t iCooper et all |2010D. Modeiling the census of ~ 
25 known Milky Way dSphs for survey incompleteness 
(|Koposov et al.l 120081 : IWalsh et al.ll2009f ). one can esti- 
mate that the number of such galaxies within the halo of 
the Milky Way could be several hundred, although with a 
signifi cant error- bar (jTollerud et al.ll2008HKoposoy et alJ 
[20091) . Nearly 50 dSph are now known around M31, 
where the whole area has been studied, but not to such 
low luminosities. Empirical constraints on the dark halo 
masses of these objects are required in order to say 
whether these discoveries represent a step toward con- 
sistency with standard cosmological models. 

The structure of the paper is as follows: we first out- 
line our specific motivation, deriving reliable velocity dis- 
persions in faint dSph galaxies. We then introduce the 
Bootes I galaxy and our observational approach, which 
is designed to test the reliability with which we can de- 
rive accurate radial velocities of faint stars. In section 4 
we describe our data reduction procedures, which are de- 
veloped to ensure use of the full information content in 
the raw data. In Section [S] we describe our enhanced 
procedure for fitting the stellar spectra and deriving pre- 
cise, accurate radial velocities, as well as our approach to 
understanding the errors. In section 6 we quantify how 
we detect radial velocity variations. Section 7 compares 
our results with available literature studies. In Section |8] 
we describe our statistical methodology to determine ve- 
locity dispersions from kinematic data, and apply this to 
the Bootes I dwarf galaxy. In Section [9] we conclude our 
study with some discussion of the astrophysical implica- 
tions. 

2. MEASURING VELOCITY DISPERSIONS IN DSPH 
GALAXIES 

Several factors complicate determinations of dSph 
galaxy masses. The most basic stems from the fact 
that we can measure only line-of-sight velocities of dSph 
stars. As a result, any estimation of a dSph's to- 
tal mass must contend with the fundamental degener- 
acy between mass and anisotropy in the velocity distri- 
bution (e.g.. [Wilkinson et al l[200l iKlevna et al.ll2002f ). 
The second basic constraint, as more generally, is that 
we very rarely see an outer declining dispersion profile 
(or rotation curve), allowing determination of a "total" 
mass. A third key limit comes from having to suppose 
the dSph can be considered in dynamic equilibrium - 
an assumption whose validity is certainly not obvious 
for those extreme systems currently within a few tens 
of kpc of the Galactic centre. Accepting the equilib- 
rium assumption, and using just radial velocities, sev- 
eral recent studies showed that the most robust mass 
it is possible to estimate for faint d warfs is the mass 
enclosed within the half light radius (iPenarrubia et all 
[20081 IWalker et al1[2009l : IWolf et al.ll2010D . However, es- 
timates of even this quantity meet with observational 
hurdles: one is that the ultra-faint dwarf galaxies have 
very small number of stars (~ 10^""*), so for some ob- 



jects the number of targets available for spectroscopic 
observation s with even the larges t telescopes is limited to 
a few tens ([Simon fc Gehal[2007[ ) of predominantly faint 
stars, for which radial velocity errors may be of the order 
of or larger than the velocity dispersion of the system. 
Another difficulty of measuring the small velocity dis- 
persions for the faint dwarf galaxies using small number 
of stars is that it still unclear how much they can be 
affec ted by the binary sta r s which must be in the sam- 
ple ("Olszewski et al.' 'l996'; 'Hargreaves. Gilmore. AnnanI 
1996; Minor ct al. 2010; McConnachic & Cote 20l3). A 
perennial limitation in all these studies is data quality 
- are systematic and random uncertainties in the data 
quantified as precisely as is claimed? 

We have two ambitions in this study. The first is to 
improve our knowledge of the Bootes I dSph galaxy. Its 
high scientific interest is introduced in the first subsec- 
tion below. Our second aim is to improve our knowl- 
edge of how reliably and accurately we can measure ra- 
dial velocities of dSph member stars. Our observations 
have been designed to allow us to derive radial veloci- 
ties with sufficient precision, and with sufficiently well- 
known accuracy, to resolve the intrinsic dispersion in a 
typical faint dSph galaxy, to properly understand our 
measurement errors, to quantify the fraction of objects 
with variable radial velocities, and hence to measure the 
intrinsic velocity dispersion with reliable precision. This 
study provides a template for how reliably, and how ac- 
curately, velocity dispersions can really be determined in 
the faintest dSph galaxies. These ultra-faint dSph have 
very few apparently bright member stars, so one is forced 
to determine precision velocities for very faint stars from 
moderate spectral resolution, moderate signal-noise ra- 
tio spectra, with the relevant velocity accuracy being at 
most a few percent of one pixel in the observed spectrum. 
The stars additionally tend to be very metal-poor, weak- 
ening absorption line strengths. 

Kinematic studies of ultrafaint dSphs are especially 
challenging due to the fact that the measured veloc- 
ity dispersions the typical errors of individual veloci- 
ties, and the the expected contribution of binary stars 
are all similar, at a few km/s. The previously studied 
SDSS dSphs h ave reported velocity dispersions as sma ll 
as ~ 3 km s"^ (ISimo n fc Geha]l20Q7l: [Martin et al.l l2007l). 
iMcConnachie fc Co te (2010) demonstrate that for such 
cold dispersions, it is difficult to disentangle contribu- 
tions from random( i.e., reflective of the underlyin grav- 
itational potential ) and binary-orbital motions. It may 
even be that some of the apparently small extreme ultra- 
faint systems are not (currently) dark matter dominated, 
but are dissolved star clusters, or tidal debris. In such a 
case the velocity dispersion could be substantially lower 
than ^ 3 km s^^. Could we detect such an absence of 
dark matter? 

In this situation it is critical to account accurately 
for velocity errors — when propagated through calcula- 
tions, optimistic/pessimistic errors result in over-/under- 
estimates of the true velocity dispersions. That is, how 
important are the uncertainties on the error bars? We 
explore the severity of this risk with some very sim- 
ple Monte Carlo simulations, in which we draw veloc- 
ity samples of 100 stars (slightly larger than our present 
Bootes I sample), in order to minimize sampling er- 
ror and thereby to isolate bias due to inadequate res- 
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Figure 1. Accuracy of simulated velocity dispersion estimates as 
a function of the ratio of measurement errors to the true velocity 
dispersion. Filled circles represent cases in which the measure- 
ment errors are known perfectly. Open circles and open squares 
represent cases in which the errors used in the analysis are un- 
derestimated (optimistic) by a factor of 0.5 and 0.75, respectively. 
Open triangles and crosses represent cases in which the adopted 
errors are overestimated (pessimistic) by a factor of 1.25 and 1.5, 
respectively. 

olution) from a Gaussian distribution with true disper- 
sion (Jtrue, and then add to each velocity a "true" error 
drawn from a second Gaussian with dispersion Stme- In- 
real observations the true error is unknown; we therefore 
measure the velocity dispersion of each artificial sam- 
ple after adopting velocity errors kStme- The factor k 
is held constant for a given sample and represents the 
accuracy of the adopted errors: fc = 1 if the adopted 
errors are accurate; A; > 1 for unduly-pessimistic er- 
rors; fc < 1 for over-optimistic errors. We repeat the 
simulation for values k = 0.5,0.75,1.0,1.25,1.5, and 
Strue/o'true = 0.1, 0.2, ...,1.5; the latter sequence lets us 
examine accuracy as the true velocity dispersion domi- 
nates, or is dominated by, the measurement errors. 

Figure [T] displays the squared ratio of the measured to 
the true velocity dispersion, (^^easuredl '^tme^ as a func- 
tion of 5true/ f^true- Plotted points represent average val- 
ues from 10^ trials at each accuracy level. We recover 
the required outcome, that if measurement errors are 
known perfectly, intrinsic dispersions can be measured 
accurately, even if they are dominated by measurement 
errors. However, as the velocity error becomes a signif- 
icant fraction of the true dispersion, deductions quickly 
become unreliable for even modestly misjudged errors. 
There is a particular danger of grossly underestimat- 
ing velocity dispersions with pessimistic errors: when 
5true ~ o'true, the measured dispersion is only 65% (20%) 
of the true value if the errors are overestimated by a fac- 
tor of fc = 1.25(fc = 1.5). 

These results emphasize the caution necessary when 
measuring and interpreting small velocity dispersions. 
Two prerequisites for obtaining reliable results are 1) 
sufficient resolution such that typical velocity errors are 



smaller than the velocity dispersion, and 2) accurate es- 
timates of the velocity errors. We designed this study 
with that lesson in mind. 

3. THE Bootes I DWARF GALAXY 

Bootes I was one of the first new dSphs discovered 
using the SDSS photometric survey, by IBelokurov et al.l 
([2006). Bootes I has a number of interesting prop- 
erties, but seems representative of the group of 
newly-discovered intrinsically-faint dSph galaxies 
which are (reasonably) far from the Galactic cen- 
tre. Many recent studies are available. Photometric 
and stellar pop ulation studies have been c ompleted 
bv IHughes. W allcrstcin & BossJ] (|2008[) : Ide Jongj 



(let al. 20081): TO kamoto (2010) ; RRLy r ae va riability 
studies by [Siegei (2006, ): DairOra et al.l (1200^ : Spec- 
troscopic abu nd ance stud ies by iMunoz et al.l (l2006l): 
iMartin et all (I2007D: iNorris. Gilmore. Wvse et al.l 
(I2008D: INorris. Yong. Gilmore et al.l (I2010D : 

INorris. Wy sc. Gilm ore et al.l (120101 ): and a HI 21cm 
search bv [Bail vn fc Ford I ( 20071) . ^Of spec ific relevance 



here, two kinematic studies are available. IMunoz et al.l 
([2006D obtained spectra of red giant branch (RGB) 
candidates over part of Bootes I, selected from SDSS 
DR4. They used the WIYN telescope and the Hy- 
dra multifiber spectrograph. Their data (for a 7 
member-star sub-sample) yielded a systemic velocity 
of 95.6 ± 3.4km/s, and a central velocity dispersion of 
6.6 ± 2.3 km/s. They derived ctq = 9.0 ± 2.2km/s from a 
12 member-star sample, a v a lue w h ich was adopt e d in a 
later analysis bv lWolf et all (|2010t ). IMartin et all (|2007D 
observed candidate Bootes I red giants from SDSS 
(DR4) with Keck/DEIMOS, finding a mean velocity 
of 99.9 ± 2.4 km/s, with central velocity dispersion 
cr = 6.5!t]^'3 km/s for their final sample of 24 member 
stars with small {6^ < 10 km/s) velocity u ncertainties. 
From their kinematics, IMunoz et al.l ()2006D deduced a 
mass of l.ltol x lO^M©. IWolf et al.l (l20Tol) adopt the 
high velocity dispersion of 9.0 zt 2.2 km/s , by a ccepting 
the superset of the data of Muf ioz et al.l (|2006l ) . From 
this they deduce a correspondingly higher mass. 

In summary, Bootes I is some 65 kpc (m — = 19.07) 
distant, has absolute luminosity My^Q = —5.9, is devoid 
of HI, has somewhat elliptical (ellipticity =0.2) morphol- 
ogy, has a half-light radius of 240pc, has an apparently 
exclusively old metal-poor stellar population, with a sig- 
nificant blue straggler sub-population, has a mean [Fe/H] 
metallicity of -2.55, has an intrinsic abundance disper- 
sion with formal a = 0.45, with a range in [Fe/H] of at 
least 1.7dex, and has at least one member star which has 
[Fe/H]= -3.7. 

Bootes I is ideal for a more detailed kinematic study. 
Observationally, Bootes I has kinematic data avail- 
able from two quite different spectrographs (one fibre- 
coupled, one slit), so that both an internal and an ex- 
ternal test of the accuracy of our data is possible. It is 
at intermediate distance from the Galactic centre (70 kpc 
Galacto-centric) , far enough that tidal effects need not be 
dominant, yet close enough that observations far enough 
down the giant branch to obtain useful statistical sam- 
ples are feasible, in spite of its (interestingly) low in- 
trinsic luminosity. It has extremely low surface density, 
so it is a clear test case for galaxy formation models, 
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none of which naturally creates very large, very low sur- 
face brightness, galaxies. It has extremely low metallicity 
stars. The abundance results suggest Bootes I is a sur- 
vivor of a true primordial system, quite likely forming 
prior to reionisation, enhancing its interest. Its velocity 
dispersion is reported as being at least 6km/s, a value 
which is feasible to measure to high accuracy. It seems 
to have an extremely high apparent mass to light ratio, 
again enhancing its intrinsic interest. It is visible from 
the VLT. The primary challenge is that it is large on 
the sky. Our present study, involving a single FLAMES 
field, as with those in the literature, covers only the cen- 
tral half-light radius of Bootes I. 

3.1. Target Selection 

For these observations, we selected target stars with 
a sufficiently wide range of apparent magnitude that a 
single observation of the brightest stars would provide 
similar signal-noise ratio to the final signal-noise ratio 
of the fainter stars, integrated over 16 identical single 
exposures. Thus, by comparing the actual repeatabil- 
ity of velocities for the brighter stars, we could verify 
if the deduced accuracy calculated for the faintest stars 
was reliable. Additionally, by summing the brighter 
star exposures, we could readily test for any system- 
atic floor in delivered accuracy at some high signal- 
noise value. We were fortunate to have access to the 
VLT FLAMES spectrograph, since that, being fibre- 
fed, reduces substantially the extra complexity of slit- 
centering errors, which unavoidably plag ues slit spec- 
troscopy (jHargreaves. Gilmore. et al]|1994D . 

We used the FLAMES spectrograph at the 8.2-m 
Kueyen (VLT/UT2) telescope at Cerro Paranal, Chile, 
to acquire spectra of individual stellar targets in Bootes I. 
Observations took place in service mode during 2009 
February-March in fulfillment of ESO Programme 182. B- 
0372A (PI: Gilmore). For this programme we use d 
FLAMES in UVES-Fibre mode (|Pasquini et al.|[200l : 
that is, 130 fibres fed the medium-resolution Giraffe 
spectrograph while eight additional fibres fed the high- 
resolution UVES spectrograph. The Giraffe spectra sam- 
ple the region SISOA- 9375A, in LR08 setup, includ- 
ing the prominent Call-triplet (CaT) absorption feature, 
with resolving power R ~ 6500. The UVES spectra re- 
sults will be described elsewehere (Monaco, Gilmore et 
al, in prep). Here we present results based primarily on 
the Giraffe spectra. 

In order to maximise the probability of observing 
Bootes I members, we selected stellar targets based 
on colors and magnitudes of the stars. The left-hand 
panel of Figure [2] displays the colour- magnitude diagram 
(CMD) of all stars (including foreground) within 12 ar- 
cmin of the center of Bootes I, from SDSS photometry 
and with ma gnitudes corrected fo r extinction using the 
dust maps of ISchlegel et al.l ()1998() . Recall that the half- 
light radius of Bootes I is greater than 12arcmin - we 
are sampling only the very inner regions in this (first) 
study. Bootes Fs low luminosity presents a challenge 
to spectroscopic studies, as the paucity of bright RGB 
candidates limits sample size. As the RGB of Bootes I 
together with a few blue horizontal branch stars(BHB) 
are clearly vis ible on the left panel o f Figure [^see also 
Figure 1 from iBelokurov et al.l (|2OO60 ) the targets were 
selected using a CMD mask covering the location of RGB 



and BHB stars. The middle panel of Figure [5] plots the 
CMD for our FLAMES targets, which include stars up 
to two magnitudes fainter than the horizontal branch. 
The UVES targets are all brighter than r ^ 18, while 
the limiting magnitude of r ~ 22 for our Giraffe targets 
represents a compromise between quantity and quality of 
the spectra in our sample. As noted above, we retain sev- 
eral (relatively) bright stars of low a priori membership 
probability as a key part of our internal quality checks. 

The second compromise involved the trade-off between 
target numbers and adequate data to test our delivered 
velocity accuracy. In order to obtain useful spectra for 
the faintest targets, and to implement our test of achiev- 
able precision, we observe only a single field with a sin- 
gle target configuration centered on Bootes I, building 
signal-to- noise (S/N) ratios by repeated science expo- 
sures. [The alternative would have been to reallocate 
those fibres allocated to brighter target stars after one or 
a few integrations, to increase observed numbers.] Our 
adopted strategy is well-suited for service mode obser- 
vations, given the excellent stability of the FLAMES in- 
strumentation. Over the six weeks of our observing pro- 
gramme, we obtained 16 individual science exposures, 
typically of 45 minutes each. The total exposure time on 
our field was 11.5hours. 

4. DATA REDUCTION 

We first summarise the sequence of our data process- 
ing, before providing a detailed description below: 

• Default basic processing using the ESO pipeline 
(bias subtraction, flatfielding, wavelength calibra- 
tion, extraction); 

• Wavelength recalibration of each extracted spec- 
trum using sky emission lines; 

• Combining sky fibres for the determination of the 
master sky spectrum for each frame and subtract- 
ing from individual extracted objects; 

• Combining the individual frames into the co-added 
spectra. Co-added spectra are used then to reject 
cosmic rays from individual spectra; 

• Spectral fitting of the co-added spectra. We deter- 
mine the best-fit template and approximate veloc- 
ity; 

• Spectral fits of individual (not co-added) spectra 
using the best-fit template. 

The initial data processing was done using the giraf- 
3.8.1 pipeline provided by ESO with some modifications 
and bug fixes described below. The important bug we 
discovered and fixed in the pipeline was related to the 
computation of the variance spectra, which had been in- 
correctly scaled by the pipeline. Another modification 
that we applied to the pipeline was related to the ex- 
tracted but not-rebinned spectra. Our goal was to min- 
imise the number of rebinning steps of the data, so we 
needed to produce non-rebinned extracted spectra as well 
as the corresponding wavelength solutions. As the orig- 
inal ESO pipeline does not provide these, we modified 
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Figure 2. Left panel: The color-magnitude diagram of all objects classified as stars with good photometry in a 12armin square area 
centred on B ootes I, from SDSS. The line shows the fiducial RGB sequence for the metal-poor globular cluster M92 (from 
IClem et al.i.20 08^ , placed at the distance of Bootes I dSph Middle panel: The color-magnitude diagram of candidate stars targeted 
for VLT spectroscopy. Right panel: The density of SDSS stars around the center of Bootes I is shown in grayscale, the targeted stars are 
shown by red points. For comparison, the half-light radius of Bootes I is 12.5arcmin. 



our version of the pipeline to output the necessary infor- 
mation, in un-rebinned pixel space. All the further dis- 
cussion of the data reduction will take into account that 
each individual spectrum is in its native pixel space, and 
so each spectrum is on a different wavelength grid. By 
contrast, the common approach has all the data rebinned 
to a common wavelength scale, and often a common dis- 
persion, prior to all subsequent anlayses. 

4.1. Wavelength calibration 

The wavelength calibration was done using the stan- 
dard Thorium- Argon arc spectra, which were taken dur- 
ing the day, not in parallel with the nighttime observa- 
tions. According to the GIRAFFE user manual this cal- 
ibration will deliver a precision limited only by spectral 
resolution and signal-noise ratio, with a delivered instru- 
mental floor at very high spectral resolution of 30m/s. 
After analysing the extracted spectra, we noticed that 
the sky lines had systematic velocity shifts of the order 
of 1— 3km/s. The left panel of Figure [3] shows these 
offsets measured using 3 isolated sky-lines (8310.7246A, 
8415.2422A and 8452.2656A) as a function of fibre ID 
for one of the GIRAFFE frames. A waving pattern of 
velocity offsets versus the fibre IDs is clearly visible. In 
order to correct for these velocity offsets in each fibre 
we used all (around 100) bright sky lines simultaneously. 
For each fibre we first subtracted the continuum from the 
sky or from the stellar spectra, then we fitted the resid- 
ual spectrum in pixel space by a sum of nxLSFs (Line 
Spread Functions) at the locations give n by the cata- 
logue of sky lines from iHanuschikl (|2003[ ) after applying 
the polynomial correction to the wavelengths. 
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Although the number of lines is high and the number 
of parameters is high also, the fit may be done easily via 
sparse matrix operations. As a result of this procedure 
we derive the best polynomial correction to the spectra 
vo,vi,V2- The results of determination of vo,vi, V2 for all 
the fibres on one particular GIRAFFE frame are shown 
in the right panel of Figure [3l Despite considerable noise 
in the measured coefficients from fibre to fibre, overall 
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Figure 3. Left panel: Velocity offsets from the ThAr arc 
wavelength fit, determined using the 3 sky lines 8310.7246A, 
8415.2422A, 8452.2656A, plotted for a single data frame, as a func- 
tion of fibre ID number Right panel: The parameters of the poly- 
nomial wavelength corrections vo,vi,V2 from Eq. [T]as a function of 
Fibre ID. Different symbols show the measurements for individual 
fibres, while the lines show the polynomial fits to these measure- 
ments. Further explanantion is in the text. 

the coefficients vary smoothly with Fibre ID (Fibre ID is 
enumerating fibers along the slit and on the CCD across 
the dispersion). In order to further reduce fibre-to fibre 
noise in the wavelength corrections we fitted the vq,vi, V2 
coefficients by a set of low-order polynomials (5th order 
for vq, 2nd order for vi and 0th order for W2). The results 
of the fit are shown as coloured lines in the right panel 
of Figure [3]) . We use those curves to determine the val- 
ues of the parameters vo,vi,V2 for each fibre, and then 
substitute those values into Eq. [T] in order to determine 
the final wavelength solution for each fibre. 

4.2. Sky subtraction 

In each Giraffe frame we allocated 12 target fibres to 
sky, with sky areas chosen to have no object detected 
by SDSS nearby. In order to obtain an average sky spec- 
trum for each frame, we first applied the corrections from 
Sec. 14.11 then rebinned the sky spectra, with twofold 
oversampling, from each individual fibre to the common 
wavelength grid. We then multiplied each of the indi- 
vidual sky spectra by a low-order(2) polynomial of wave- 
length to ensure that all rebinned sky spectra had the 
same photon-count normalisation with wavelength, and 
median-combined the rebinned sky spectra, creating an 
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'average' sky spectrum. For the error spectruni of the 
combined sky spectrum we used the median absolute de- 
viation (MAD), sc aled by 1.4826^/|"/sqr t(number of sky 
spectra) , following iRousseeuw fc Croux I ()1993[ ) . 

The average sky was then subtracted from each indi- 
vidual fibre on the frame. In order to perform this sub- 
traction we had to rebin the average sky spectrum back 
onto the original pixel wavelength grid of each fibre. Af- 
ter subtracting the scaled sky spectrum from each fibre 
we also added the variance of the mean sky spectrum to 
the variance spectrum of each fibre to take into account 
the uncertainty of the sky determination. 

4.3. Combining spectra 

Since each object in our sample was observed from 9 
to 18 times, we need to co-add the individual spectra. 
Although we did not use the co-added spectra for the 
radial velocity determinations, we needed the co-added 
spectra for certain specific tasks, for example, cosmic ray 
rejection in individual exposures. 

In order to combine repeated science spectra for a given 
star we first made a zeropoint wavelength correction to 
correct for the varying radial velocity component due to 
the Earth's motion, then rebinned each individual sky- 
subtracted spectrum to a common wavelength grid, as 
previously done for raw sky spectra. Then we followed a 
procedure similar to that which we used in determination 
of the average sky spectrum. Since our spectra were not 
flux corrected, before combining the individual spectra 
we used 2nd order polynomials to multiply all individ- 
ual spectra, to ensure that each had the same continuum 
level across the wavelength range. Then we median com- 
bined these spectra. For variance spectra we used, as 
before, the median absolute deviation(MAD), scaled by 
1.4826-\/^/sqrt(number of observations). 

Having the co-added spectra for each object allows us 
to identify bad pixels/ cosmic rays/outliers in each in- 
dividual (not co-added) spectrum. In order to do that, 
we interpolated the co- added spectra to the wavelength 
grid of each individual observation, and then we masked 
out those pixels which were more than 4 sigma above or 
6 sigma below the median spectrum. Strictly speaking, 
that procedure may introduce biases if the spectrum is 
highly variable from exposure to exposure - in that cases 
the variable lines may be masked out, but we visually 
checked all the spectra and did not see any inappropri- 
ately masked spectral features. 



l2007tlKoleva et al.|[2009l) . Methods based on direct pixel 
fitting provide more realistic error-bars and give a better 
way to treat multiple templates and continuum levels. 

In this section we briefly describe the pixel-fitting 
method we use, which i s similar to the ones de scribed 
in Rix fc Whitd (ll992l):jCappellari fc EmsellemI (|2004) : 
iChilingarian et al.l H20071 ): iKoleva et all (|2009[ ) but with 
a few differences. 

5.1. Stellar library 

An important ingredient for direct pixel-fitting meth- 
ods is the library of template spectra. We decided to use 
the lib rary of synthetic spectra provided bv lMunari et al.l 
(|2005f ). The spectra in the library cover a large range of 
stellar parameters -2.5 <[Fe/H]< 0.5, [a/Fe] = 0,0.4, 
3000 < Teff/l/sT < 80000, 1.5 < \og{g) < 5, < 
Vrot/^km/ s < 150. The highest resolution spectra avail- 
able from the Munari library are R ~20000, which is 
higher than the resolution of our data (R~6500), so the 
templates can be easily downsampled to our resolution. 
While synthetic spectral libraries have obvious limita- 
tions and we do not expect a perfect match to the ob- 
served stars, the range of stellar parameters covered in 
the library is large and is significantly better than what 
we can achieve with stellar libraries of high resolution 
spect ra for rea l stars such as EL ODIE (jPrugniel et al.l 
I2007D : see also iKirbv et al.l (|2008D . This is particularly 
true given the extremely low metallicities of the B ootes I 
member stars (|Norris. Wvse. Gilmore et aLll2010[ ). 

5.2. Fitting synthetic templates to real spectra 

Having the stellar template library we construct the 
model of each observed spectrum as 

Model{\ i, V, {pj}) = P(A) • T, (a (1 + ^)) (2) 

where Ti (A) is the i-th template spectrum from the tem- 
plate grid convolved with the appropriate Line Spread 
Function LSF{X), v is the radial velocity of the object 

and P{X) — Pj ■ A-' is the normalising polynomial of 

3=0 

low degree, (N-1), which takes into account the lack of 
flux calibration of our spectra, as well as any imprecision 
in our location of the continuum of the template spectra. 
Having the spectral model one can then compute the 
value by summing scaled residuals over pixels. 



5. VELOCITY DETERMINATION AND SPECTRAL 
FITTING 

The standard appr oach used for the measurement 
of stellar velocit ie s (iHargreaves. Gilmore. et al.l 119941: 
Koch et all [20071: ISimon fc Gehal 120071: IWalker et al.l 



20071: iBaumgardt et al.l 120091: iGeha et al.l 120091: 
Leaman et al. l l2009l:IKoch et al. Illoof has been to cross- 



correlate against a template spectrum (jSimkinl 119741 : 
iTonry fc DavisI 119791 ) . Although simple and computa- 
tionally fas t, cross-correlation is known not to perform 
optim ally ()Rix fc White! 119921 : iCappellari fc EmsellemI 
l2004f ). In fact direct pixel- fitting methods have been 
widely employed for more than a decade in spectroscopic 
studi e s of unresolved stellar populations (Rix fc White! 
119921: iCappcUari fc Emsellem ,2004; .Chilingarian ct al.i 



( Spk - Model{Xk,i, v, {pj}) 



k ^ 



ESpk 



(3) 

where Spk and ESpk are the observed spectra and vari- 
ance spectra respectively, and the Xk are the wavelengths 
of the pixels of the extracted spectra. Model{Xk, ■■■) is 
the evaluation of the synthetic stellar model spectrum at 
the wavelengths of the individual pixels. It is important 
that the original observed spectrum is neither rebinned 
nor interpolated in any way, as this would lead to corre- 
lated noise and information loss. 

Equation [3] defines the or log-likelihood of our 
model. In order to find the estimated radial velocity 
and the best-matched template we need to minimise the 
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Figure 4. Examples of fitting synthetic model spectra to co-added spectra for several stars of different magnitudes (from top to bot- 
tom): SDSS J135922.59-I-143300.7 (r = 18.51), SDSS J140001. 53+142154.2 (r = 19.35), SDSS J140002.44-I-142249.1 (r = 20.62), SDSS 
J135951.07-f 143049.8 (r = 21.52). Black lines show the observed co-added spectra, while the red lines show the best fit model spectra. 



with respect to all relevant paranieters. In fact we 
are not interested in values of the coefficients of the nor- 
malizing polynomials P{X), so we can marginalise over 
them and derive the joint probability distribution of just 
the template identification i, and the radial velocity, v. 
With a simple analytical computation we can derive that 
probability distribution: 



Prob{i,v) = {dct{M)) 2 exp 



MY 



(4) 



where Y is the vector having length N, such that : 



ESpk 



(5) 



and M is a symmetric square NxN matrix, such that 



ESpi 



(6) 



Equation 2] describes the joint probability distribu- 
tion of an identified template and an associated target 
star radial velocity Prob{i,v). In practice we evaluate 
this probability for a grid of plausible radial velocities 
\v\ < 500km/s (this is effectively a uniform prior on ra- 
dial velocity) and for our grid of templates (also assuming 
a uniform prior). This 2-dimensional probability distri- 
bution can then be marginalised over template or veloc- 
ity. The marginalisation over template can be used to 
determine the maximum a-posteriori (MAP) estimate of 
the velocity, (argmax^, Prob{v)), and of the velocity er- 
ror. The other marginalisation identifies the best fitting 



template spectrum, and therefore an estimate of stellar 
parameters. 

There are a few important points to appreciate about 
the fitting procedure used here. First, the data do not 
have to be rebinned to either linear or logarithmic wave- 
length scales; therefore all the information as well as the 
noise properties are preserved. Second, we do not need 
to perform continuum subtraction from either the object 
or the template, which is advantageous since that is al- 
ways a poorly defined task; instead we rely on the contin- 
uum shape from the synthetic spectrum, wi t h add itional 
polynomial modifications (cf iKoleva et al.l ()2009( ) for a 
discussion of continuum fitting) . 

Figure |4] shows the result of the spectral fitting proce- 
dure applied to the co-added spectra of four of the stars 
in the Bootes I sample. Figure [5] shows the resulting 
velocity probability distribution for the same four stars. 
The velocity precision ranges from 0.25 km/s to 7km/s. 
For the brighter stars the probability distribution of the 
velocity is approximately Gaussian, while for the fainter 
star it is clearly asymmetric. For some of the faintest 
stars in our sample the probability distribution is even 
multi-modal. 

5.3. Application to Bootes I stellar spectra 

For each star, as a first step we apply our fitting 
procedure to the co-added spectra from each exposure. 
The purpose of using the co-added spectra here is pri- 
marily to determine the best-fit template for each ob- 
ject and the subsequent velocity estimate. Although 
this fitting is motivated by optimizing the radial veloc- 
ity accuracy, and not primarily to measure the stellar 
astrophysical parameters T^s, \og{g) and [Fe/H] from 
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Figure 5. Un-normalized velocity posterior probability distribu- 
tions derived from the co-added spectra, for the same four stars 
as shown in Fig. (4] It is clear that, especially at faint magni- 
tudes/lower signal- noise, the velocity probability distribution is not 
Gaussian. 

our template fit, the values of these parameters deter- 
mined from the best-fit templates are reasonable, as il- 
lustrated by Figure IHl This Figure shows how the mea- 
sured parameters correlate with radial velocities and col- 
ors of the stars. The left panel of Figure [S] shows the 2d 
distribution of radial velocities and metallicities of the 
stars. From the plot we see that the velocity peak at 
100 km/s related to the Bootes I dwarf galaxy is pro- 
duced by stars with [Fe/H] < —1.5, which is what we ex- 
pect from the metallicity of Bootes I ()Martin et al.|[2007l : 
iNorris. Gilmore. Wvse et aLll2008[ ). In the middle panel 
of Figure |6] we show the 2D distribution of radial veloc- 
ities and log(g). We see that the contamination (from 
foreground Milky Way stars) typically has high surface 
gravity log{g) > 4.0 - exactly what we expect from disk 
dwarfs. The right panel of the figure shows the correla- 
tion between determined effective temperature and (g— r) 
colour. Although of relatively low significance a correla- 
tion still can clearly be seen. Overall we conclude that 
the parameters from the template fit are reasonably well 
determined using our fitting procedure and the adopted 
spectral library. 

As the next step we use the best fit templates as well 
as the initial velocity estimates for the final fitting of 
each individual (i.e., not co-added) observation of a given 
object in the same way as described above. The velocity 
estimate from the fit to the co-added spectra, plus or 
minus 50 km/s, is used as a uniform prior on the radial 
velocity. That is, we make an assumption that the radial 
velocity of a given object in a single observation does 
not vary more than 50 km/s from the value measured 



from the co-added spectrum. From this step we obtain 
probability distributions of the velocities for each epoch 
observation of each object. We can use these repeated 
velocity measurements for several purposes: first in order 
to assess our measurement errors, and second to check for 
possible binarity /variability of the radial velocities (see 
Section [5]). We end with a total sample of 112 stars for 
which we have derived radial velocities. 

5.4. Checking the Derived Uncertainties 

As described in the previous section, from our spectral 
fitting procedure applied to either individual or co-added 
spectra we determine the probability distribution of the 
radial velocity of a given star Prob{v), the MAP (maxi- 
mum a-posteriori) estimate of the velocity and the error 
of that velocity. One of the most important checks on 
the validity of our results — as well as a check on the ef- 
fectiveness of our reduction method — is the confirmation 
that the velocity errors that we determine from individ- 
ual observations of specific (non-variable) stars are not 
systematically larger or smaller than the scatter between 
individual repeated measurements. That is, have we met 
our ambition of deriving correct and reliable uncertain- 
ties? 

Figure [7| compares the standard deviation determined 
from repeated measurements to the mean error deter- 
mined by our fitting procedure. In this plot we expect to 
see a one-to-one correlation, except for three cases: one, 
where there is intrinsic variability in the radial veloc- 
ity; two, where there is some remaining systematic error 
which varies from observation to observation; or three, 
in the case where measurement errors from individual 
exposures are significantly larger than the one-pixel res- 
olution of the spectrograph ('^15 km/s). In this latter 
case undersampling causes the probability distribution 
of the velocities to become significantly non-Gaussian, 
so the mean fitting error is an inadequate representation 
of possible variation of the velocities. 

With reassuringly few exceptions, the data points in 
Figure [7| scatter about the one-to-one line. There is a 
slight apparent systematic tendency for the data points 
to be scattered more above the one-to-one line than be- 
low it, expected given that the error-bars on the standard 
deviation are significantly asymmetric. The one-to-one 
relation extends in general down to a precision of a few 
hundred m/s, confirming that our derived error-bars are 
reliable down to 0.1— 0.2 km/s. In the next section we 
show that several points lie above the one-to-one line not 
because we have underestimated our errors, but rather 
due to intrinsic variability of their radial velocities. 

Since we were somewhat, albeit pleasantly, surprised 
at obtaining velocity precision substantially better than 
Ikm/s, as confirmed by repeated measurements, we 
checked whether this may be some artifact of the al- 
gorithm due to the sky lines. For example, if a sys- 
tematically poorly subtracted sky line affecting a spe- 
cific wavelength caused velocities to cluster near some 
fixed value, this situation could mimic small variances in 
repeat observations. A piece of circumstantial evidence 
against this hypothesis is that the nearly one-to-one rela- 
tion of the error-bars versus the standard deviation from 
repeated measurements (Fig. [7]) contains a set of stars 
with very different radial velocities, while it seems un- 
likely that sky-line artifacts are present at all possible 
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Figure 6. The stellar parameters determined from our synthetic spectrum template fit. The left panel shows the 2D histogram of radial 
velocities and metallicities of the best-fit templates for all observed stars. It is clear that the velocity peak at ^ lOOkm/s, the systemic 
velocity of Bootes I, occurs at [Fe/H] < — 2dex, which is ~ the metallicity of Bootes I, while the background (halo, thin and thick disk) has 
[Fe/H] > —1.5. The middle panel shows that most Bootes I stars have low surface gravities, as expected for giants, while the contaminating 
stars from the Milky Way foreground have almost exclusively log{g) > 4.0, as expected for dwarfs. The right panel shows the correlation 
of the g-r colour vs the effective temperature of the best-fit template. This correlation demonstrates that the effective temperatures are 
reasonably well determined. 

highest precision(< Ikm/s from individual exposures), 
there is a possibhty that the systematic errors due to 
template mismatch (e.g., stellar spectra looking system- 
atically different from spectra in our grid) are dominating 
the error budget. 

6. VARIABILITY ANALYSIS 

The goal of this section is to use our velocities, and 
velocity uncertainties, to estimate which of the observed 
stars show significant intrinsic variability, so that we may 
remove those stars from the analysis of the velocity distri- 
bution in Bootes I. We do this by testing the plausibility 
of two hypotheses for each object, the first hypothesis 
being that the object does not show significant varia- 
tion in radial velocities, and the second hypothesis being 
that the object does show evidence for velocity variability 
with amplitude greater than Ikm/s. 

For each star and each of N observations of each star 
our spectroscopic fitting procedure provides us with the 
probability distribution of the radial velocity Probi{v) 
associated with that observation. In general these prob- 
ability distributions are not always Gaussian, especially 
in the low signal to noise cases, so in the following analy- 
sis we will avoid making any assumptions of Gaussianity. 

For each star we assume that there is a certain intrinsic 
distribution of radial velocities "^{v — Vsys), where Vsys 
is the systemic velocity of the star. If the radial velocity 
of the source is not varying then \E'(w — Vsys) is a delta 
function. The likelihood of the observations, assuming 
that they randomly sample ^(w), can then be written as 
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Figure 7. Comparison of the mean velocity error determined from 
the fitting procedure and the estimated standard deviation of the 
repeatedly measured velocities. Error bars show 68% confidence 
limits for the estimated standard deviations. The solid black line 
shows the one-to-one relation. The grey line shows the approximate 
velocity where the one-to-one relation between the mean error and 
the estimated standard deviation is expected to fail due to non- 
Gaussianity of the probability distributions. 

radial velocities. The most compelling evidence against 
this hypothesis is that, since we have radial velocity mea- 
surements spread over a month of observations, the vari- 
ation of barycentric corrections is roughly 3km/s, which 
is much larger than our claimed precision and scatter be- 
tween repeated measurements; this would not be possible 
if our measurements were driven by some earth-velocity 
sky-features or artifacts. 

As a final cross-check of our radial velocity error-bars, 
we compared them with the theoretical estimates of min- 
imum possible velocit y uncertainty, based on Fisher- 
matr ix like arguments ()Murphv et al.|[2008l : iGriest et alj 
l201Cft. This comparison suggests that our error-bars are 
not unrealistic, and are typically very close to that min- 
imum velocity uncertainty, being above it only for a few 
of the faintest stars. This confirms that our velocity pre- 
cision is not unreasonable. 

A last remark is that , while the random errors seem to 
be correctly determined for the brightest stars with the 



L=l[ hfiv 



Probi{v) dv. 



(7) 



It is clear for example that if all the individual Proht {v) 
are Gaussians, and if the object's velocity is assumed to 
be constant, (ie, "^{v — Vsys) is a delta function), then the 
likelihood L will be a Gaussian centered at the weighted 
mean of the centers of the Gaussians {Probi(v)), i.e. in- 
dividual velocity measurements. 
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Figure 8. The probability of detected radial velocity variability 
for all 112 observed stars, computed using the Bayes factors for our 
two hypotheses: HI, there is no intrinsic variability in radial ve- 
locities; and H2, there is intrinsic variability with amplitude larger 
than Ikm/s. The grey line shows the cut in Bayes factor, which 
is placed such that stars above the line favour the H2 hypothesis, 
and are intrinsically variable. There are 11 such stars. The blue 
diamond indicates a known RR Lyrae star. 
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Figure 9. Radial velocity variation of the star SDSS 
J135951. 33+143905. 8 , identified as an RR Lyrae variable by 
IDall'Ora etall lf200^) and Siegel (2006). The left panel shows the 
observed radial velocities as a function of the Julian date, while 
the right panel shows the o bserva t ions p hased with the photomet- 
ric period P=0.3119d from lSiegell i20m) . 

However, we do not make these Gaussian assumptions, 
and can adopt a more general approach. 

In order to assess the velocity variabihty ol a given 
object we evaluate the likelihoods of our two hypotheses: 
HI, that the velocity of the object does not change, that 



is '^{v 



^sys 



) = 6iv 



5); and H2, that the velocity 



is varying and the velocities are distributed uniformly 



between Vsys 
the scatter: 



s and Vs 



■ s. where s is the measure of 



^ if If — V 



if \v 



sys 



< S 

> s 



(8) 



For each hypothesis H1,H2 we can compute its likeli- 
hood using Eq. [71 and marginalise over the parameters 
Vsys and s, adopting uniform priors. The ratio of the 
likelihoods of two hypotheses marginalised over the pa- 
rameters becomes a Bayes factor, B - which is a measure 
of the relative plausibility of the two hypo theses peffrevsl 
[TOBlt IKass fc R,aftervl[T995l : iTrottal [200I . We calculate 
this for all 112 stars in our final sample. 

Figure |8] shows the probabilities of variability calcu- 
lated for all 112 of our stars, derived from Bayes factors, 
as a function of stellar apparent magnitude. The grey 
line shows the location of the probability threshold (see 
[Kass fc Raftcrv,,1995i) which we use to identify 11 stars 
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Figure 10. Radial velocity variations for two stars in our sample, 
one indentifified as variable, one not. Left panel: The radial veloc- 
ity observations as a function of time for the fairly bright star in 
our sample, J135940.18-f 142428.0. Note the vertical scale. This is 
shown to illustrate the accuracy with which we are determining ra- 
dial velocities, right panel: Radial velocity variation of star SDSS 
J140003.32-I-142851.4. 

with significant evidence for radial velocity variation. 

The blue diamond point, located significantly above 
the grey line, corresponds to one of our target stars. Af- 
ter (re-)discovering its variability, we realised this star is 
a known RR Lyrae variable, from publishe d photometric 
studies (|DairOra et al.l[200l ISiegell 120061) . This recov- 
ery does confirm that we can indeed correctly identify 
stars with variable radial velocity. Figure [9] shows the 
observed (as a function of observation time), and phased 
using the photometric period radial velocity variations 
for this RR Lyrae star. 

Figure [TU] illustrates the observed velocity data as a 
function of time for one bright, high signal- noise ratio 
star not detected to be velocity variable, and for another 
star with variability Bayes factors above our variability 
detection threshold. 

In short, the discussion above, and the evidence in the 
shown in these figures, demonstrate that we are able to 
calculate a Bayes factor probability that each observed 
star is consistent with being radial velocity non- variable 
over the time in which we have observations. While pass- 
ing this test says little about much longer period velocity 
variability, it does allow us reliably to identify stars for 
which we have radial velocity data, but whose kinemat- 
ics cannot be included in a dynamical analysis at face 
value. Applying that criterion, we restrict further anal- 
ysis to 100 stars, with adequate velocity measurements 
and error bars, and no evidence for velocity variability. 

7. COMPARISON WITH PUBLISHED KINEMATIC STUDIES 
OF Bootes I 

Two observational kinematic studi es of Bootes I 
(Mufioz et al. 2006; Martin et al' ll2007D. and two other 
analyses (jW olf ct al. 2010; Norris. Wyse. Gilmore et al.l 
l2010t iNorris. Yong . Gilmore et al..i201C)l ). have been pub- 
lished, providing information which we may compare 
w ith our result s. 

iMunoz et all (|2006D used the WIYN telescope and 
the Hydra multifiber spectrograph to measure radial 
velocities for 58 candidate member stars, all brighter 
than 5=19. Considering the strength of the Mgb ab- 
sorption features, they classify fully 30 of these as gi- 
ants. Further considering position, velocity and line- 
strength, they identify seven stars located within lOar- 
cmin of the centre of Bootes I to define the mean veloc- 
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ity and velocity dispersion of Bootes I, finding a mean 
velocity 95.6 ± 3.4km/s, and central velocity dispersion 
6.6 ± 2.3km/s. With these values they further identify 
a total of 12 stars as "potential Str members", 11 being 
within a radius of ISarcmin of the centre of Bootes I, the 
12th at 27arcmin distance (the half-light radius is 12.5ar- 
cmin). These 12 candidate members provide a higher 
mean velocity, 98.4 ± 2.9km/s, and a higher dispersion 
of (To = 9.0 ± 2.2km/s; this higher dispersion, as the au- 
thors note "corroborating the apparent increase of the 
velocity dispersion with radius" . The authors also note 
they have two further "likely high- velocity members" . 
The full set of 14 stars provides a velocity dispersion of 
(Jo = 14.6 ± 3.0 km/s, suggesting "a possibly dramatic in- 
crease of the Boo velocity dispersion with radius" . The 
12-star sample, with its Bootes I vel ocity disper s ion o f 
(To = 9.0±2.2km/s, is that adopted bv lWolf et aD (2010) 
in their dynamical analysis. We compare our results for 
9 stars in common in the left panel of Figure [TT] Our 
results in this study a re inconsisten t with the velocity 
dispersions reported bv lMufioz et al.l ()2006[) . and we rule 
out their suggested rapid radial increase in velocity dis- 
p ersion (see be l ow). 

iMartin et al.l (|2007l ) observed 96 candidate Bootes I 
red giants from SDSS (DR4) with Keck/DEIMOS. They 
identified a sample of 24 member stars, each with a ra- 
dial velocity deter mined with an accuracy smaller than 
6 km/s (cf. Marti n et al.l I2007L for details). From this 
sample they find a mean velocity of 99.9 ±2.4 km/s, with 
central velocity dispersion a = &.bt^'\]iTii/s for their fi- 
nal sample of 24 stars with small velocity uncertainties. 
We observed, by design, 27 stars in common with the 
Martin et al. sample. 

We compare our measurements for the 27 stars in 
common with the radial velocity measurements by 
IMartin et al.l (|2007| ) in the rightmost two panels of Fig- 
ure [TTJ The zero-points of our radial velocity measure- 
ments are consistent. The mean velocities, for the com- 
mon 27-star member sample, are Vkeck = 99.6±1.7km/s, 
and VvLT,27 = 101.2 ± 2.0km/s, in excellent agree- 
ment. The mean velocity from our full sample is Vvlt = 
101.8±0.7km/s. The right hand panel of Figure [TT] com- 
pares our der i yed v elocity errors with those reported by 
IMartin et al.l ()2007[ ). It is apparent our derived single- 
velocity errors are a factor of about two smaller than 
those of IMartin et al.l (|2007[ ). at high signal- noise ratios 
in both studies. We may also do a very crude check on the 
accuracy of the quoted v elocity uncertainties . Compar- 
ing stars in common with 'Marti n et aP (|2007[ ). the veloc- 
ity difference, in units normalised b y the quadrat u re sum 
of our quoted errors and those of (Martin et~all (|2007D . 
has a sigma of 2.7. This is robust evidence the combined 
errors are underestimated. The results of Figure [7] sug- 
gest we have correctly calculated our uncertainties. In 
addition, as we discuss in the section below, our derived 
velocity dispersion for Bootes I is significantly below that 
published by other studies. We interpret this difference, 
in essence a difference between fibre-fed and slit spec- 
trographs, as reflecting the inherent precision limits of 
velocity determination using multi-slit spectrographs. 

For completeness, we note the studies 

of ^Bootes I by iNorris. Wvse. Gilmore et al.l 

(|2010f) and INorris. Yong. Gilmore et al.l ()2010[ ). 



INorris. Wvse. Gilmore et all ()2010D obtained spec- 
tra of candidate RGB members within a 1 degree radius 
of Bootes I, using the AAT-|-AAOmega facility. Their 
study was designed to identify candidate members of 
Bootes I at large distances from the centre, for analysis 
of the chemical abundance distribution, and to allow 
subsequent more detailed study for kinematics, so 
they had, by design, relatively low velocity precision 
(10 km/s). Nonetheless, they identify candidate mem- 
bers out to 60arcmin (5 half-lig ht radii, cf their Fig 6) 
from t he centre of Bootes I. Norr is. Wvse. Gilmore ct a P 
( 2010j ) observed 5 stars in common with IMartin et al.l 
1 2007f ). and report agreement in mean velocity within 
3 km/s and a dispersion in velocity differences for those 
5 stars of 2.3 km/s (in spite of their suggested 10 km/s 
accuracy). They further report work in preparation 
confirming one of their candidates (Boo-980) to be 
an extremely metal-poor giant at 3.9 half-light radii 
from the centre of Bo5tes I, and with radial v elocity 
99.0 ± 0.5 km/s. INorris. Yong. Gilmore et al.l ()2010( ) 
report on a follow- up high dispersion UVES stu(i y of on e 
star identified by INorris. Wyse. Gilmore et al.l (|2010( ). 
This star, Boo-1137, is also extremely metal-poor. 
Of relevance here, it lies 24arcmin (2 half-light radii) 
from the centre of Bootes I, and has radial velocity 
Vuves — 99.1 ± 0.5 km/s. Both these velocities are in 
remarkable agreement with the mean velocity of the 
Bootes I system, considering their large Bootes-centric 
distance, and hint at a low velocity dispersion in the 
outer very metal-poor populations in Bootes I. 

8. THE VELOCITY DISPERSION(S) IN Bootes I 

Equation [7] gives us a way to determine the best ve- 
locity estimate from several velocity measurements. Af- 
ter removing those stars which we suspect to be vari- 
able in radial velocity, we model the remaining 100 star 
sample of velocities assuming that each velocity is con- 
stant (i.e. '^{v — Vsys) is a delta- function) . Thus we 
derive the probability distribution of the systemic veloc- 
ity for each star by multiplying radial velocity probabil- 
ity distributions from separate measurements Prob{v) — 
Y[ Probi(v). Those probability distributions are 

observations 

in most cases very close to Gaussians, so in the following 
analysis we assume that we have the velocity estimate 
and the Gaussian error-bar for each star. Table [T] lists 
those velocities, and other information, for each star. 

Figure [T2l shows our resulting velocity distributions for 
two different subsets of stars. The top panel shows the 
distribution of velocities for the 37 stars which are likely 
members, ie those with log (5) < 3.5, [Fe/H]< —1.5, not 
showing any variability according to our Bayes factor 
criterion and with small errors of their radial velocity 
measurement < 2.5 km/s. The bottom panel shows 
the corresponding velocity distribution for all 100 non- 
velocity variable stars in our sample. The velocity dis- 
tributions illustrate two important points. First, the ve- 
locity peak due to Bootes I member stars is quite evi- 
dent. Second, it is clear that the velocity dispersion of 
the stars in Bootes I is significantly smaller than 6 km/s, 
or 9 km/s, the glob a l dispersions measured previously 
iMufioz et ail (|2006t ): IMartin et ahl £2007). In the top 
panel especially, the bulk of the distribution looks similar 
to a Gaussian with ~ 3 km/s dispersion, while in the bot- 
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Figure 11. Left panel: Comparison of our velocity measurements for 9 stars in common with the velocity measurements of lMufioz et al.l 
||200 6:) . The line shows the expected one-to-one relation. The data sets are apparently not consistent withing the published errors, and 
suggest a velocity-dependent scale error. Middle panel: Comparison of our velocity measurements for 27 stars in common with the velocity 
measurements of Martin ct al. (2007). The line shows the expected one-to-one relation. The red points identify stars which may have 
variability in radial velocities (see Section |6]l. Right panel: Comparison of our radial velocity errors with those of .Martin et al., (,2007 ) . The 
thick black line shows the one-to-one relation, while the grey line shows the y = 3x relation. 
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Figure 12. The distribution of stellar velocities in Bootes I. 
The black line shows the distribution of velocities estimated using 
the Epanechnikov kernell lEp ancchnikov"1969l: IWand fc JoneslllOQSI ) 
with a bandwidth of 1.5 km/s, the grey line shows a standard his- 
togram with bin size of 1.5 km/s. The red and blue lines are over- 
plotted Gaussians with sigma of 3 and 6 km/s; respectively; 6 km/s 
is the smaller of the previously published determinations of the 
Bootes I internal velocity dispersion. The top panel shows the ve- 
locity distribution for our sample of 37 stars which are highly prob- 
able Bootes I members, i.e., those with [Fc/H]<-1.5, log{g) <3.5, 
small velocity error ay <2.5km/s and no significant velocity vari- 
ability. The bottom panel shows the velocity distribution for all 
100 of our stars with non- varying radial velocity 



torn panel the less restrictively selected sample shows, in 
addition to the low dispersion 'core', rather pronounced 
higher velocity tails. 

In order to assess the distribution of radial velocities 
we fit the observed velocity distribution with two differ- 
ent models: one where the velocity distribution in the 
galaxy is represented by a single Gaussian, and the sec- 
ond with the velocity distribution being the sum of two 
Gaussians with the same mean. To perform these fits, 
we follow the standard Bayesian approach. We write 
down the probability distribution as a function of the 
template metallicity, template \og{g) and radial velocity 
as the mixture of distributions for the background (with 
fraction f^g and for Bootes I: 

Fe Fe 

Fe 

(1 - hg) Pbooi-^] Pboo{l0g{g)) Pboo{v){9) 




Figure 13. Top panel: The probability distribution of the internal 
Bootes I Gaussian velocity dispersion, determined from an IVICMC 
fit to our velocity data for our full 100-star non-variable sample, 
when the velocity distribution is assumed to be consistent with a 
single Gaussian. Middle panel: the MCMC fit Gaussian, with dis- 
persion 4.6 km/s, overlaid on the kinematic data. Lower panel: to 
illustrate that the derived MCMC fit is robust to data selection, we 
show the derived Gaussian wi th d ispersion 4.6 km/s overlaid on the 
subset of 37 stars from figure 1121 those which are highly probable 
Bootes I members, ie those with [Fe/H]<-1.5, log(g) <3.5, small 
velocity error ay <2.5km/s and no significant velocity variability 

This technique minimizes the subjectivity involved in our 
selection, for illustrative purposes, of the 37-star subsam- 
ple noted above, where we subjectively imposed an as- 
trophysical prior, and which we restricted to the small 
subset of highest precision data. 

We assume that the probability distributions ff)g[§], 
Ptg{log{g)), Pbg{v), Pbooi^], Pbg{hg{g)) are Gaussians 
with different means and standard deviations. For the 
radial velocity distribution of stars in the dwarf Pboo{v) 
we assume that it is either a single Gaussian or the 
sum of two Gaussians with different dispersions but the 
same mean. Having the probability distribution de- 
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Figure 14. Top two panels: The probability distributions of the 
internal Bootes I Gaussian velocity dispersion, determined from 
an MCMC fit to our 100-star non-variable velocity data, v^hen the 
velocity distribution is assumed to be consistent with two Gaus- 
sians. The top panel shows the probability distribution for the 
lower dispersion component, while the bottom panel shows the 
probability distribution for the higher dispersion component. The 
MCMC analysis allocated 70% of the stars to the 2.4km/s dis- 
persion component, and 30% of the stars to the 9 km/s dispersion 
component. Lower panel: The corresponding two-Gaussian distri- 
bution overlaid on the kinematic data. 

fined in Eq. ^ we perform a standard Markov Chain 
Monte Carlo (MCMC) sampling of the available param- 
eter space in order to determine the posterior probabil- 
ity distribution for the parameters of Pboo{v). We use 
the standard Metropolis-Hastings algorithm (Metropolis 
et al. 1 953; Hastings 19 70) imple mented in t he pymc 
package ()Patil et al.|[20lol ): see, e.g.. lNea]| ()1993[ ) for a re- 
view of the MCMC method. 

The posterior probability distribution for the velocity 
dispersion of Bootes I for the single Gaussian hypothesis 
is shown in Figure 1131 The velocity dispersion estimate 
is then 4.6^g gkm/s. This Gaussian dispersion is shown 
overlaid on the kinematic data in the middle panel. To 
illustrate graphically that our results are robust against 
sample selection, in the bottom panel we show the same 
Gaussian distribution overlaid on the subset of 37 high- 
probability members with excellent data. The fit is ac- 
ceptable, but far from excellent. 

If we make the assumption that the Bootes I veloc- 
ity distribution consists of two Gaussians [in addition 
to the fit to the background], where one Gaussian has 
higher dispersion than the other, the posterior prob- 
ability distributions for the velocity dispersions of the 
two components are those shown in Figure 1141 The ve- 
locity dispersion of the lower dispersion component is 
then 2.4lQ gkm/s, while the velocity dispersion of the 
other component is not very well determined, but is 
around 9 km/s. The fraction of stars belonging to the 
higher dispersion component according to the MCMC 
fit is around 30%. The corresponding Gaussian distri- 
butions are overlaid on the kinematic data in the bot- 
tom panel. The mean velocity from our full sample is 



VvLT = 101.8 ± 0.7km/s. In order to assess the prob- 
ability that the Bootes I velocity distribution is indeed 
described by two Gaussians instead of one, we measure 
the likelihood ratio for these two hypotheses, which gives 
—2log{Li/ L2) = 8.06. This ratio corresponds to a ^98% 
confidence of rejecting the single Gaussian hypothesis. 

Our data also provide direct limits on kinematic gra- 
dients in Bootes I. We fit the MCMC modelhng of the 
kinematics allowing for a linear gradient in either RA 
(essentially the minor axis) or DEC (essentially the ma- 
jor axis. Our formal limits on rotation are: minor axis, 
—4 ± 9km/s/deg; major axis ± 8km/s/deg, recalling 
that our data cover a radial range of 0.2deg. We may also 
limit any radial change in the dispersion. This is more 
difficult, since the minority higher-dispersion component 
becomes dominated by sampling noise if the sample is 
subdivided by radius. Hence we fit a single Gaussian 
velocity dispersion model to the inner half (stars within 
0.06deg of the centre of Bootes I) and the outer half of 
the sample. The difference in derived single- Gaussian 
dispersion is then (outer minus inner)=— 2.1 ± 1.3 km/s. 
This is not a significant result, but hints that the disper- 
sion may be decreasing with radius. That is, the higher 
velocity dispersion component may (not significantly) be 
somewhat more centrally concentrated than is the lower 
velocity dispersion component. 

9. DISCUSSION AND CONCLUSIONS 

In this paper we explain how we have developed and 
implemented a thorough analysis of low-moderate spec- 
tral resolution [R = 6500) VLT/FLAMES steUar spec- 
tra taken in the CaT wavelength region near 860nm. We 
optimised a data reduction methodology which delivers 
very accurate radial velocities, and very reliable uncer- 
tainties on those radial velocities. We set up an opti- 
mised observational proof of methodology, targeting faint 
candidate RGB stars in the very metal-poor Bootes I 
dwarf spheroidal galaxy. By making 16 individual re- 
peat observations over six weeks, with a wide dynamic 
range in each observational data set, we have reached 
several goals: most importantly, we have been able to 
properly assess the errors of our radial velocity measure- 
ments; our delivered radial velocity precision, for faint 
extremely weak-lined stars, is better than 1 km/s for the 
combined exposures. Second, we have been able to iden- 
tify and reject stars that show significant radial velocity 
variability. 

Comparing our derived velocity dispersi on, and indi- 
vidual velocities, with li terature studies (|Munoz et al.l 
I2006t iMartin et alll2007[ ). shows that previous studies 
have substantially overestimated the velocity dispersion 
of Bootes I. It is possible that earlier studies underesti- 
mate their velocity errors, and hence overestimate the 
velocity distribution which is deconvolved from those 
errors. Stable fibre-fed spectrographs, including espe- 
cially VLT-I-FLAMES, when complemented with an ap- 
propriate observational strategy, and suitably sophisti- 
cated data processing, are able to deliver precise, reliable 
and accurate radial velocities with sufficient precision to 
resolve the intrinsically very low velocity dispersions ev- 
ident in the low luminosity dSph satallite galaxies. By 
exploiting the spectrograph stability to build integrations 
on times from days to years, we can detect many radial 
velocity variables, whose unrecognised presence would in- 
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flate erroneously a derived velocity dispersion. Impor- 
tantly, from repeated observations, we can quantify re- 
liably our velocity accuracy. We are currently applying 
this observational technique to several other dSph galax- 
ies. 

We have useful radial velocity measurements for 100 
non-variable RGB candidate stars, all within one half- 
light radius (12.5arcmin, 240pc) of the centre of Bootes I. 
Approximately 60-70 stars are likely members from our 
full sample. Implementing a general MCMC analy- 
sis, which includes separate Bootes I-member and back- 
ground distribution functions of our derived stellar pa- 
rameters log(g), [Fe/H], and of the member and back- 
ground velocity distributions, we show that the distri- 
bution function of stellar radial velocities in Bootes I 
which we measure can be described in two ways. The 
less likely is that the distribution is Gaussian, with a ve- 
locity dispersion of 4.6lQ g km/s. The more likely is that 
the distribution consists of two components: a "colder" 
component, containing 70% of the member stars, which 
has a projected radial velocity dispersion of 2.4^05 km/s, 
and a "hotter" component, containing 30% of the mem- 
ber stars, which has a projected radial velocity dispersion 
of 9 km/s. The data favor, with 98% confidence, the two 
component model. 

Our data also provide direct limits on kinematic gradi- 
ents in Bootes I. Our formal limits on rotation are: minor 
axis, — 4±9km/s/deg; major axis 0±8km/s/deg. Simi- 
larly, we may limit any radial change in dispersion. This 
is more complex, as fitting the full 2-component with a 
radial term becomes sample-size limited. However, fit- 
ting a single Gaussian dispersion to the inner half-radius 
and outer half-radius provides a formal, statistically in- 
significant, radial decrease in dispersion of 2.1±1.3km/s. 
That is, there is a non-significant hint that the appar- 
ently higher dispersion component is more centrally con- 
centrated than is the whole sample. 

We consider now the possibility that the higher veloc- 
ity dispersion component in Bootes I kinematics is an 
artifact of unresolved binaries with variable radial veloc- 
ities. This paper is based on the data taken during a 
period of one month. While it definitely allows us to re- 
ject stars with significantly variable radial velocities on 
hour to week timescales, we do not of course identify 
all binaries. A hypothetical substantial population of 
binaries with velocity amplitudes of order 10 km/s, and 
periods much longer than a month, might produce spu- 
rious wings in the radial velocity distribution which we 
would identify as a hotter compon ent. F igure 1111 limits 
the plausibility of this speculation: [Martin ct al. (2007) 
observed Bootes I in 2006, we observed Bootes I in 2009 
- comparison of their velocity data with ours for 27 stars 
in common identifies perhaps three stars with velocity 
near that of Bootes I, and with sufficient velocity vari- 
ability to populate the distribution function wings. This 
is inconsistent with the results of our MCMC analysis, 
namely that some 30% of stars populate the apparently 
higher velocity dispersion component. Thus taking the 
data of lMartin et al.l (j2007[ ) at the face value, binary vari- 
ability is unlikely to be the cause of the higher-dispersion 
component we detect in our sample. 

What can we deduce about the mass of Bootes I? There 
has been recently a flurry of studies on means to deter- 



mine some useful mass parameter to represent a dSph 
galaxy, where limited radial velocity data are available, 
concentrated in the central regions. We noted above that 
the most robust mass that can be estimated for faint 
dwarfs is the mass enclosed within the half light radius 
(for example, among very many studies. Walker et al.l 
lloollWolf etalll20ia) . These several methods m essence 
determine M1/2 = where is a character- 

istic mass inside the galaxy's half-light radius R1/2, 
is the (Gaussian) dispersion of the line of sight radial ve- 
locities, and /3 is a factor of value 2.5/G to 3/G, with 
G the Newtonian constant. For a given object, clearly 
this mass scales as the square of the velocity dispersion. 
Equally clearly, it assumes there is a single well-defined 
Gaussian dispersion. 

For Bootes I we seem to have two dispersions, with no 
robust way to associate a scale length with either. If we 
just assume that the dominant low dispersion component 
is associated with the measured half-light radius, then 
we can (very approximately) deduce an associated mass 
within that radius of 240pc. It is perhaps more useful 
to consider the r ange of dete rminations of the half-light 
mass of Bootes I. iWolf et all (|201Q1 adopt a velocity dis- 
persion for Bootes I of 9.0 km/s, which they correspond 
to a mass Afi/2 = 2.36.10'^Mq, and M/L = 1700. Fixing 
the geomet ric parameter s , but adopting the dispersion 
derived bv iMartin et al.l (|2007[ ). = 6.5km/s, corre- 
sponds to a mass lower by a factor of 1.9. Adopting our 
low dispersi on of (jy = 2.4 km /s provides a mass lower 
than that of IWolf et all ()2010 D by a factor of 14, and a 
corresponding M/L = 120. The range of these numbers, 
depending entirely on the quality of the kinematic data, 
clearly illustrates both the need for excellent quality data 
to allow useful study of the faint dSphs, and the consider- 
able caution which should be applied to extant analyses 
of the masses of very low luminosity dSph galaxies. 

9.1. Two populations in Bootes I: some speculation 

Our data suggests that the Bootes I dwarf spheroidal 
has both cold and hot stellar components. Although the 
data do not completely rule out a single stellar compo- 
nent with velocity dispersion of 4.6 km/s, that scenario 
is not favored by our data. In order to understand the 
origin of the two components we checked whether the 
components differ in radial distribution or metallicity, 
but we were not able to find any significant differences. 
This is not a strong constraint, given i) that our data 
extend over only one half-light radius from the centre of 
Bootes I, and ii) our metallicities are unable to resolve 
the extremely low abundanc es of Bootes I members. 
Norris, Wyse. Gilmore ct al.! (|2010D show (their fig 17) 
that the mean abundance of Bootes I is [Fe/H]= —2.5, 
with a tail down below [Fe/II]= —3.5. These abun- 
dances are below the bottom range of our template spec- 
tra, and beyond our ability to resolve in this study. At 
face value two velocity components correspond to two 
scale lengths. Might Bootes I have an extended higher 
dispersion component? The available Subaru photome- 
try (Okamoto, 2010) shows no indication of an extended 
'envelope' structure, but is not a strong constraint be- 
yond 2 half-light radii, as the stellar surface density is 
extremely low. We did note in Section |7| that radial ve- 
locities are available for two extremely metal-poor mem- 
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bers of Bootes I, at distances of 2.0 and 3.9 half-light radii 
from the centre. Both velocities are within 0.5kni/s of 
the systemic velocity. While only two stars, they are 
consistent with an extremely extended, very metal-poor, 
low velocity dispersion component in Bootes I. These two 
stars were however selected for detailed analysis because 
they were suspected of being very metal-poor. One may 
not deduce that the more metal-poor stars in Bootes I 
are the more extended. These two stars do however pro- 
vide some evidence against significant rotation or tidal 
warping of Bootes I, even out to 1 kpc. 

Perhaps more interesting, and consistent with a veloc- 
ity distribution function which is not a single compo- 
nent, is the radial distribution of the 16 member stars 
wi th [Fe/H] abundances derived by, aii d listed in Table 3 
of. INorris. Wvse. Gilmore"etaII mWi ). If we divide the 
16 stars into equal inner and outer groups of eight stars 
based on radius from the centre of Bootes I, the inner 
and outer groups have mean abundances [Fe/H]=— 2.30± 
0.12, and [Fe/H]=-2.78±0.17 respectively These differ 
at the 2Aa level. 

An alternative interpretation is to question the valid- 
ity of fitting Gaussians to th e radial velocity distribution 
function. In simple models ()Gerhardlll993D strong devi- 
ations from a Gaussian distribution of radial velocities 
may arise from the velocity anisotropy of the stellar pop- 
ulation. Given our lack of understanding of how dSph 
galaxies form, there is no physical basis for assuming a 
single isotropic velocity distribution. There are two good 
reasons to doubt the single gaussian assumption for the 
very faintest dSph. Firstly, just how does one populate 
an extremely low-density distribution with scale length 
250pc, out to apparently 4 scale lengths, lOOOpc, with a 
velocity dispersion below 4km/s? Dispersing a central 
star formation region in a very shallow (cored?) dark 
matter potential may be feasible. Or it may not. Such a 
process will inevitably generate very radially biased or- 
bits. One alternative speculation is to merge several star 
forming regions, none of which need have been centrally 
located, during first formation of the dark matter po- 
tential which we call Bootes I. This might well generate 
highly tangentially biassed velocity distributions, where 
one might anticipate the more tangentially-biased orbits 
to be more centrally concentrated. Both processes might 
happen, so that the velocity distribution may well be bi- 
or multi-modal in anisotropy, with significant radial gra- 
dients in this anisotropy. The least likely expectation is 
that the velocity distribution function is isotropic - there 
is no obvious physical process which could generate such 
a distribution function at such low densities and such 
large scales. 

A plausibility argument of relevance to this speculation 
is that the lowest luminosity dSphs are systematically 
less round than are more luminous galaxies - although 
Bootes I itself has ellipticity 0.2, so is relatively round. 
This flatness may be generated by an isotropic kinematic 
distribution in a flattened potential, or equally by an 
anisotropic kinematic distribution in a spherical poten- 
tial. It is often assumed that dSph shapes correspond 
to the shapes of the dominant dark matter distributions. 
Since one of the very few things we suspect we know 
about the low luminosity dSph is that mass does not 
follow light (more correctly, light does not trace mass, 
kinematics do), the shape of the luminosity distribution 



may reflect kinematic anistropy, not mass anisotropy. In 
that case radially anisotropic distributions will appear in 
projected radial velocities as a distribution which is more 
"cuspy" than is a Gaussian, while tangential anistropy 
will appear more extended than is a Gaussian ([Gerhard! 
Il993t IBinnev fc Tremguni [2008). The combination of 
these effects can look rather like our measured Bootes I 
velocity distribution in Figure [Ml 

Thus, one may interpret our statistically preferred 
two-Gaussian radial velocity distribution function in two 
ways. One option is that there is a two-component 
structure in Bootes I, with future observations at large 
Bootes-centric distances required to identify the char- 
acteristics of the higher-velocity dispersion population, 
which is a minority in the inner regions, and which must 
have a hugely extended radial scale length. Alternatively, 
Bootes I has a velocity distribution function which re- 
flects its formation, and which is a combination of a ma- 
jority population with very radially-biased orbits, and a 
minority population with very tangentially biased orbits. 
Both are represented by a single, measured, radial scale 
length. The more radially-biased population contains the 
most metal-poor stars. 

If this speculation is valid, one anticipates that fu- 
ture precision determination of the velocity distribution 
function in the most flattened dSph galaxies (Hercules, 
e = 0.6 and UMa II e = 0.54) will find neither is consis- 
tent with being a single component Gaussian, and will 
find a radially-variable shape for the distribution func- 
tion of radial velocities. 

Most importantly, and our conclusion, it is evident 
that radial velocity data of very high precision, and of 
extremely well-defined accuracy, are necessary to make 
progress in defining the kinematics and masses of the 
faintest dSph galaxies. 
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4750 -0.5 4.5 

5250 -1.5 4.0 

5000 -2.5 2.5 

4000 -2.5 3.5 

5000 -1.5 4.5 

4750 -1.5 1.5 

5000 -2.5 1.0 

5250 -2.0 2.5 B 

4500 -0.5 4.5 

XXX 

XXX 

5000 -2.5 1.5 B 

6000 -1.0 4.5 

5000 -2.5 1.0 B 

4750 -2.5 1.0 B 

XXX 

5250 -2.0 2.5 B 

5000 -2.0 2.5 B 

5000 -2.5 2.0 

5000 -2.5 4.5 

5250 -1.5 3.0 

5000 -2.5 2.0 B 

5000 -2.0 1.5 B 

4500 -0.5 4.0 

5000 -2.0 1.5 B 

5500 -2.0 1.0 B 

5000 -2.0 2.0 

5250 -1.0 4.5 

5000 -2.0 2.5 B 

5000 -2.0 1.5 

4750 -1.5 1.5 

4000 -2.5 1.0 

7000 -2.0 2.5 

4500 -2.5 1.0 B 

5000 -1.5 4.0 

6000 -2.5 4.5 

5000 -2.0 1.0 B 

5000 -2.5 3.5 

5000 -2.5 1.5 B 

5000 -1.0 2.5 

4000 -2.0 1.0 

5000 -2.5 4.5 

8000 -2.0 3.5 

5000 -2.0 1.5 B 

5000 -1.5 2.5 

4750 0.0 4.5 

5250 -1.0 4.5 

4750 -1.0 4.5 

5250 -2.0 3.5 

5000 -2.5 1.0 

7000 -2.5 4.5 

5000 -2.0 1.0 B 

9000 -1.5 3.5 

6000 -1.5 3.0 

5500 -2.0 4.0 

4750 -2.5 1.0 B 

5000 -1.5 4.5 

4750 -1.5 2.0 

5000 -2.5 2.5 

5000 -2.0 3.0 B 

5000 -2.5 1.0 B 

4750 -2.5 1.5 

6000 -2.0 3.0 

5000 -2.5 2.0 B 

5000 -2.5 1.0 B 

5000 -2.5 1.5 B 

5000 -2.5 1.0 B 

5000 -2.5 1.0 

5000 -2.0 2.0 B 

4500 -2.5 1.0 B 
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Table 2 
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r 


RV 


RV error 




Tefl 


[Fe/H] 


log(5) 


Bestflag 






[dcg] 


[dcg] 


[mag] 


[mag] 


[km/s] 


[km/s] 




[K] 
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TT Ar\r\-\ 1 tr o i i /lotrcc r\ 


21U.U4oU 
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21.0 


114.41 


2.07 


0.62 


5000 


-2.5 
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B 


/ i 




01 n nc: 1 n 
2IU.U0IU 


1/1 AQC\A 

14.4oy4 


21.3 


20.9 
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2.61 


0.75 
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-2.0 


4.5 
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J i4UUiz. 41 + 14,5.127.0 


m n n cr 1 "7 

21U.U517 


14.5576 


20.6 


20.3 


-90.00 


1.40 


0.52 
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-2.5 


4.5 




/ O 


J14UUlz.yz+i4ooil. i 


1 n n Q 
2IU.UO00 


14.5533 


19.6 


19.0 
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0.31 


0.53 
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-2.0 
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/4 
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1 n OiT^Kd 
2IU.UOOO 


1/1 AOQ A 

14.4oo4 


18.7 


18.0 


44.92 


0.16 


0.52 


4750 


-2.0 


4.5 
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J14UU14.Di + 14oyoU. i 


1 n Old 1 1 
2IU.U0II 


1/1 av^^QK 
14.D0O0 


21.3 


21.0 


99.91 
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0.55 
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-1.5 


4.0 
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01 n cxa A n 
21U.Uo4U 
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14.oo42 


18.2 
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0. 12 


0.31 
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-1.0 


3.5 
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0.84 
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-2.5 
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J14(J(J1d.1oH-14o14d.9 


01 n n^JTo 

210.0673 


14.5297 


21.2 


20.7 


X 


X 
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X 


X 


X 
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J 140016.59+143530.0 


210.0691 


14.5917 


19.6 


19.7 


99.39 


1.21 


0.64 


7500 


-1.0 


3.0 




81 


J 14(JU1d. 02+142925. 4 


210.0693 


1 A A nn A 

14.4904 


21.3 


20.9 


97.21 


2.09 


0.55 


5000 


-2.5 


3.5 
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Ji4UU2U.o8+14d7o4.U 


m n no cr 

21U. 1)858 


14.6261 


22.1 


21.7 


102.01 


5.57 


0.90 


5000 


-2.5 


1.5 
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T-| /innoi nn i i /iQnoQ o 
J 14UU21.UU+14oy2o.2 


01 n novc; 
21U.U8 (d 


14. DODO 


19.5 


18.9 


-30.74 


0.32 


0.61 


4750 


-1.0 


1.5 
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n /innoi 0/i i i /IOcjcq q 
J 14UU21.o4+ 142000.0 


01 n nni n 
21U.UylU 


14.4315 


21.0 


20.5 


90.41 


1.38 


0.59 
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1.5 
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J 14UU22.1U+1400O0.2 
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14.D4oy 


19.5 
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0.55 
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01 n nno cr 

210.0935 


14.5575 


19.6 


19.2 


98.47 


0.44 


0.49 
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-2.0 


1.5 
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87 


Ji4UU2o.oo+142DU7.y 


21U.Uy72 


14.4356 


21.9 


21.3 


92.79 


2.30 


0.56 
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1.0 
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0.84 
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0.53 
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1.0 


rS 
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01 n 1 c\(\A 
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14.57D2 


18.8 


18.2 


on TA 
29.74 


0.17 


U.5U 


5000 


1 n 
-l.U 


A C 

4.5 




92 


J 14U(J2d. 52+142919.0 


01 nil nc 

21U.11U5 


1 A A oon 

14.4809 


21.5 


21.2 


-117.49 


2.59 


0.60 


5000 


-2.0 


1.5 




93 


J 14UU2D.57+14294o.y 


01 n 11 nT 

2lU.llu7 


14.49d9 


21.2 


20.7 


108.10 


1.52 


0.53 


5000 


-2.5 


3.0 


B 


CiA 

94 


Ti/inno^; 0*7 1 i/i/ion/i 
J 14UU2D.O i -\- L44iZ\j4l.2 


01 nil on 
21U. 112U 


1 A 'rm 
14. / U12 


21.9 


21.5 


-27.44 


9.14 


1.00 


4000 


-2.5 


3.0 




95 


TI /innoT r\ A i i /looon o 

J 140027. 04+ 143830. 3 


01 n 1 1 T 

21U.1127 


14.6418 


19.6 


19.7 


105.37 


1.81 


0.55 


8000 


-2.0 


3.0 


B 


96 


TI /innoT 00 1 1 /looin cr 

Ji4UU27.2o+14o21y.5 


01 n 1 1 T7 

21U.llo7 


1/1 cr 

14.5o8o 


20.3 


19.8 


105.68 


0.79 


0.51 


4750 


-2.5 


1.0 


B 


97 
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14.5533 


22.2 
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J 14UU2o.oy+142oo2.D 


01 n 1 1 Q 
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14.09 /9 


21.4 


21.1 


98.85 


2.94 


0.65 


5000 


-2.5 


2.0 
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TI /innoo H'O 1 1 /I 1 /1 T 

Ji4UU28.7o+14oi42.7 


1 n 1 1 nT 

210.1197 


14.5285 


21.8 


21.3 


2.25 


2.25 


0.55 


5000 


-1.5 


4.5 




100 


TI /innoo no 1 1 /ioc;no 
J i4UU2o.yo+142oU2.2 


1 n 1 oni^ 
21U. 12Uu 


1/1 /II TQ 

14.41 / 


21.8 


21.4 


103.40 


3.65 


0.58 
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-2.5 


2.0 
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J14UUol.oo+14o/ lo.D 
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2IU.I0U5 
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21.9 


21.5 


-43.94 


3.64 


0.60 
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J 140031.78+142015.5 


01 n 1 00 A 

21U.la24 


14.3376 


19.2 


18.7 


6.44 


0.41 


0.49 


5500 


-1.0 


4.5 




104 


TI dnnoi ni i i AAif\o o 

J14(JUol.91+1441Uo.o 


01 n 1 oon 

21(J.lo29 


14.6856 


22.4 


22.1 


136.32 


6.38 


0.75 


4000 


-2.5 


3.5 




105 


TI /innoo 1 /I 1 1 oifOT a 

J 14(JUo2. 14+14od27.d 


01 n 1 oon 

210.1339 


14.6077 


21.9 


21.4 


164.05 


3.63 


0.56 


5000 


-2.5 


3.0 




106 


TI /innoo oT 1 1 ,<o'70c o 

J 140032.37+142735.8 


210.1349 


14.4bUU 


22.2 


21.9 


X 


X 


X 


X 


X 


X 




107 


J140032. 54+142400. 6 


210.1356 


14.4002 


19.1 


18.5 


-60.56 


0.19 


0.76 


5250 


-0.5 


4.5 




108 


J140033.07+142959.7 


210.1378 


14.4999 


19.7 


19.2 


96.86 


0.44 


0.52 


5000 


-2.5 


2.0 


B 


109 


J140033.75+142514.3 


210.1406 


14.4206 


18.6 


17.9 


-40.08 


0.13 


0.39 


4750 


-1.0 


4.5 




110 


J140035.97+142854.5 


210.1499 


14.4818 


21.6 


21.1 


-123.72 


2.82 


0.56 


5000 


-2.5 


4.5 




111 


J140037.38+142858.0 


210.1558 


14.4828 


17.4 


16.4 


106.95 


0.16 


1.00 


4500 


-2.0 


1.0 




112 


J140039.56+142827.8 


210.1648 


14.4744 


19.6 


19.7 


95.13 


3.71 


0.71 


8000 


-2.5 


3.5 




113 


J140040.90+143208.1 


210.1704 


14.5356 


21.9 


21.7 


106.83 


4.28 


0.56 


6000 


-2.5 


4.5 




114 


J140047.55+142412.0 


210.1981 


14.4033 


21.0 


20.5 


117.74 


1.62 


0.61 


5000 


-2.5 


2.5 


B 


115 


J140047.57+142630.2 


210.1982 


14.4417 


21.3 


20.9 


99.13 


2.14 


0.57 


5000 


-2.5 


1.5 


B 


116 


J140053. 19+142705.3 


210.2217 


14.4515 


16.6 


16.1 


42.30 


0.07 


0.99 


4500 


0.0 


4.5 




117 


J140057.91+142854.1 


210.2413 


14.4817 


21.7 


21.2 


103.03 


3.83 


0.57 


5000 


-2.0 


2.0 





